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DNA polymerase (DNAP) is a dual-purpose enzyme that plays two opposite roles in two different 
situations during DNA replication. It plays its normal role as a polymerase catalyzing the elongation 
of a new DNA molecule by adding a monomer. However, it can switch to the role of an exonu- 
clease and shorten the same DNA by cleavage of the last incorporated monomer from the nascent 
DNA. Just as misincorporated nucleotides can escape exonuclease causing replication error, correct 
nucleotide may get sacrificed unnecessarily by erroneous cleavage. The interplay of polymerase and 
exonuclease activities of a DNAP is explored here by developing a minimal stochastic kinetic model 
of DNA replication. Exact analytical expressions are derived for a few key statistical distributions; 
these characterize the temporal patterns in the mechanical stepping and the chemical (cleavage) 
reaction. The Michaelis-Menten-like analytical expression derived for the average rates of these two 
processes not only demonstrate the effects of their coupling, but are also utilized to measure the 
extent of replication error and erroneous cleavage. 

PACS numbers: 87.16. Ac 89.20.-a 



I. INTRODUCTION 

DNA polymerase (DNAP) replicates a DNA molecule; 
the sequence of the nucleotides, the monomer ic subunit 
of DNA, on the product of polymerization is dictated 
by that on the corresponding template DNA through 
the Watson- Crick complimentary base-paring rule 
DNAP moves step-by-step along the template strand uti- 
lizing chemical energy input and, therefore, these are also 
regarded as a molecular motor [l, Q . 

An unique feature of DNAP is that it is a dual-purpose 
enzyme that plays two opposite roles in two different cir- 
cumstances during DNA replication. It plays its normal 
role as a polymerase catalyzing the elongation of a new 
DNA molecule. However, upon committing an error by 
the misincorporation of a wrong nucleotide, it switches 
its role to that of a exonuclease catalyzing the shortening 
of the nascent DNA by cleavage of the misincorporated 
nucleotide at the growing tip of the elongating DNA [4] . 
The two distinct sites on the DNAP where, respectively, 
polymerization and cleavage are catalyzed, are separated 
by 3-4 nm [5]. The nascent DNA is transferred back to 
the site of polymerization after cleaving the incorrect nu- 
cleotide from its growing tip. The elongation and cleav- 
age reactions are thus coupled by the transfer of the DNA 
between the sites of polymerase and exonuclease activity 
of the DNAP. However, the physical mechanism of this 
transfer is not well understood 

In this paper we develop a minimal kinetic model of 
DNA replication (more precisely, that of the "leading 
strand" which can proceed continuously) that captures 
the coupled polymerase and exonuclease activities of a 
DNAP within the same theoretical framework. From 
this model, we derive the exact analytical expressions 
for (i) the dwell time distribution (DTD) of the DNAP 
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at the successive nucleotides on the template DNA, and 
(ii) the distribution of the turnover times (DTT) of the 
exonuclease (i.e., the time intervals between the suc- 
cessive events of cleavage of misincorporated nucleotide 
from the nascent DNA). The mean of these two distri- 
butions characterize the average rates of elongation and 
cleavage, repectively; we show that both can be written 
as Michaelis-Menten-like expressions for enzymatic reac- 
tions which reveals the effect of coupling explicitly. 

In our model, the kinetic pathways available to the cor- 
rect and incorrect nucleotides are the same. However, it 
is the ratio of the rate constants that makes a pathway 
more favorable to one species than to the other. Similar 
assumption was made by Galas and Branscomb [7] in one 
of the earliest models of replication. Therefore, in spite 
of the elaborate quality control system, some misincor- 
porated nucleotides can escape cleavage; such replication 
error in the final product is usually about 1 in 10^ nu- 
cleotides. Moreover, occasionally a correct nucleotide is 
erroneously cleaved unneccessarily; such "futile" cycles 
slow down replication [8]. 

We define quantitative measures of these two types of 
error and derive their exact analytical expressions from 
our model for "wild type" DNAP. Using special cases of 
these analytical expressions, we also examine the effects 
of two different mutations of the DNAP jEf- (i) "exo- 
deficient" mutant that is incapable of exonuclease activ- 
ity, and (ii) "transfer-deficient" mutant on which the rate 
of transfer to the exonuclease site is drastically reduced. 



II. MODEL 

Almost all DNAP share a common "right-hand-like" 
structure. Binding of the correct dNTP substrate trig- 
gers closing of the "hand" which is required for the forma- 
tion of the diester bond between the recruited nucleotide 
monomer and the elongating DNA molecule. The kinetic 
scheme of our stochastic model of replication is shown in 
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FIG. 1: The kinetic model of DNA replication. The chemo- 
mechanical states of a single DNAP and the allowed transi- 
tions are shown along with the corresponding transition prob- 
abilities per unit time (the rate constants). The corresponding 
rate constants for the incorrect nucleotide are denoted by the 
symbol Q (see the text for details). 



fig ([I]). The rate constants for the correct and incorrect 
nucleotides are denoted by u and 1], respectively; the 
same subscript is used in both the cases for the same 
transition. 

Let us begin with the situation where the DNAP is 
ready to begin its next elongation cycle; these mechano- 
chemical state is labelled by the integer index 1. In prin- 
ciple, the transition 1^2 consists of two steps: binding 
of the dNTP substrate and the formation of the diester 
bond. The overall rate of this step is ujf for a correct 
substrate and for an incorrect substrate. 

Occasionally, because of the random fluctuation of the 
"hand" between the "open" and "closed" conformations, 
the dNTP may escape even before the formation of the 
diester bond; this takes place with the rate constant ujr- 
If the recruited dNTP is incorrect, the hand remains 
"open" most of the time and the rate constant for the 
rejection of the dNTP is Qr (^r >> ^r)- Note that 
dNTP selection through 1^2 involves a discrimination 
between the correct and incorrect dNTP substrate on 
the basis of free energy gained by complementary base- 
pairing with the template. 

The transition 2(z + l) l{i -\- 1) corresponds to the 
relaxation of the freshly incorporated nucleotide to a con- 
formation that allows the DNAP to be ready for the next 
cycle. The rate constants for this step are Uh and fth 
respectively, for correctly and incorrectly incorporated 
nucleotides. Alternatively, while in the state 2{i + 1), 
the DNAP can transfer the growing DNA molecule to 
its exo nuclease site; this transfer takes place at a rate 
ujpf i^pf) if the selected nucleotide is correct (incorrect). 
Since Qh « and Vtpf » ujpf, the misincorporated 
nucleotide most often gets transferred to the exonuclease 
site whereas relaxation, rather than transfer, is the most 
probable pathway when the incorporated nucleotide is 



FIG. 2: "Purely" polymerizing reaction (left panel) and 
"pure" exonuclease reaction (right panel). 



correct. 

The actual cleavage of the diester bond that severs 
the nucleotide at the growing tip of the DNA is repre- 
sented by the transition 3(i + 1) l(i). For a correct 
nucleotide, ujpr » oJe indicating that the DNA is likely 
to be transferred back to the polymerase site without 
the unneccesary cleavage of the correct nucleotide. In 
contrast, for an incorrect nucleotide, Qpr << Qe which 
makes error correction a highly probable event. More- 
over, Qpr << UJprj >> ^e- 

Since the trimmed DNA is transferred to the poly- 
merase site extremely rapidly [4], each of the rate con- 
stants uje and Qe incorporate both the trimming and 
transfer. 

Interestingly, the full kinetic scheme in figH] can be 
viewed as a coupling of a purely polymerase-catalyzed 
reaction (shown in the left panel of figEJ and a purely 
exonuclease-catalyzed reaction (shown in the right panel 
of fig 12]); the transition with the rate oOpr couples these 
two reactions. 

Strictly speaking, for an "exo-deficient" DNAP [5|, 
cjg = = although the rates of forward and reverse 
transfer between the sites of polymerase and exonuclease 
activities may not be necessarily negligible. Similarly, ei- 
ther ujpf = = r^p/, or ujpr = = Qpr (or, both) can be 
the cause of "transfer-deficiency" of the DNAP. 



III. RESULTS AND DISCUSSION 

A. Distribution of Dwell time 

The DTD considered here arises from intrinsic stochas- 
ticity and not caused by any sequence inhomogeneity of 
the mRNA template [9]. For a molecular motor that is 
allowed to step backward as well as forward, we use pos- 
itive (+) and negative ( — ) signs to represent the forward 

and backward steps, respectively. For example, i/j-i {t) 

is the conditional DTD (cDTD) [10] when a forward step 

is followed by a backward step and p-\ is the probability 

of such a transition. Therefore, the DTD can be written 
as 

V^(t) =p+_v^+_(t)+p++v^++(t)+p_+v^_+(t)+p__v^__(t) 

(1) 
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Since in our model two consecutive backward steps are 

forbidden implies that p--ip {t) = 0. We calculate the 

cDTD following the standard method [11] that has been 
used successfully earlier for the calculation of cDTD for 
some other motors (see, for example, ref.[10]). 

Let P^{j,t) be the probability of finding the DNAP 
in the ja-th {/a = 1,2,3) chemical state at the j-th site 
(i.e., at the discrete position Xj {j = oo, ... — 1, 0, 1, ...oo). 
Then master equations for (j , t) are 



dt 



dP2{j.t) 
dt 



-UJfPi{j,t)^UJrP2{j + l,t)^UJhP2{j.t) 



UJePsU^l.t) 



(2) 



UJprPsU.t) 



(3) 



dt 



= UJpfP2{j, t) - {uJe + UJpr)P3{j^ t) (4) 



In terms of the Fourier transform 

CO 



iqxj 



(5) 



of P^{xj^t)^ the master equations can be written as a 
matrix equation 



jP{q,t)=M{ci)P{q,t) 



(6) 



where P(g', t) is a column vector whose 3 components are 
P^{q,t),P2{q,t),Ps{q,t) and 



M(q) 



UJpf -{uJe^UJpr) 



(7) 

with p-\-{q) = e and p-{q) = e*^*^; d being the step 
size, ie., Xj+i — Xj = d. 

Taking Laplace transform of (|6j) with respect to time 



Jo 



(8) 



the solution of the master equation in the Fourier-Laplace 
space is 



where 



P{q,s) = R{q,s)-'P{0) 



R(g,s) = si - M{q) 



(9) 



(10) 



and P(0) is the column vector corresponding to the initial 
probabilities. 



Now we define 



P{q,s) = Y,Pii<l,s) 



(11) 



which can be calculated from 



(12) 



where Coj^i are the cofactors of the R(q,s). 

The determinant of the matrix R(q,s) is a third order 
polynomial of s and can be expressed as 

|R(^, s)\ =s^^ as^ + f3{q)s + j{q) (13) 

Note that a is independent of q whereas /3 and 7 are 
the function of q. For the explicit form ([7|) of M the 
coefficients a, P{q) and j{q) are given below. 

a = CJe + Cc;/ + CJ/i + UJpf + LJpr + UJr (14) 

P can be expressed as 

m = m+P+{i-p+{q))+P+-{i-p+{q)p-iq)) (is) 

where 



/3+ UJfUJh 



(16) 



(17) 



and 



= UJfUJr. (18) 

Similarly, 

l{q) = 7+(l - P+{q)) + 7+-(l - P+{q)p-{q)) (19) 
where 

7+ = UJeUJfUJh + LdfUJh^pr (20) 

and 

7^ = UJeUJfUJpf + UJeUJfUJr + UJf UJprUJr (21) 

For convenience, we define the 2x2 diagonal matrix 



P(q) = 



p+iq) 

p_(g) 



the column vector 



*(s) = - 

s 



1 -P++V^++(5) -P+-V^+-(5) 

1 -p_ + V^_+(s) -p__V^__(5) 



(22) 



(23) 
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and the 2x2 matrix 



Equation (|27|) can be re expressed as 



p- 



-is) 
-is) 



(24) 



where tp±±{s) are the Laplace transforms of the cDTDs 
P{q, s) and cDTD are related by the equation 



P{q,s)=pUl-i^{s)p{q))-'^is) 



(25) 



where pg is the vector of initial conditions. For exam- 
ple, = (10) corresponds to the given condition that 
the motor has taken the initial step in the forward (+) 
direction. 

Thus, in principle, if one can calculate P{q^ s), one can 
use the relation (|25]) to solve for and, then taking 

inverse Laplace transform, obtained '0±±(t). To calculate 
P(g, s), one has to use an appropriate set of initial con- 
ditions consistent with the defnition of the dwell times. 
The set Pi = 0,P2 = l^Ps = ensures that first step 
is taken forward. In other words, in our calculation, we 
start the clock by setting it to t = when the DNAP 
reaches the state 2 at j from state 1 at j — 1. Therefore, 

in this case p = = p ^. Corresponding to this initial 

condition, we now define, 



P+{q,s) = P{q,s] 

{Pi(0)=0,P2(0) = l,P3(0)=0} 

and, from equation (|25]) . we get [11] 

I - p+{q)p++i^++{s) 



(26) 



sP^{q,s) 



(27) 



1 



sP^{q,s) 



ao H-a+p+((7) 



(28) 



{p-iQ)=0} 



where 



1 



ao 



Hence, 



1 -p++^++(<5) -p+-V^+-(<5) 



1 -p++^++(<5) -p+-^+-(5) 



ao 



(29) 



(30) 



and 



p+_v;+_(s) 



Therefore, next we obtain 



ao + <2+ — 1 



ao 



1 



(31) 



directly 



{p-(q)=o} 



sP^{q,s) 

from (p!2|) and, by comparing it with equation (|28|) . find 
out the expressions for ao and a+; substituting these ex- 
pressions for ao and a+ into equations (|3Q|) and (|3T]) we 

get p++?/;++(s) and p+_?/;+_(5), respectively. 

Using the same initial condition, from equation (p!2]) . 
we get 



~ + S{a - UJrjl - P-{q)) + /3(0) - (1 - P-{q)){uJeUJpf ^ Ue^r + LOpr^r) 

(32) 

Therefore, 

1 



sP+{q,s) 



^ + + s{i3{0) + /3+ + /3+_) + 7+ + 7+_ - (g/3+ + 7+)p+(^) ^33^ 



Comparing equation (|33]) with equation (|28|) we identify ao and a+ and substituting these expressions for ao and a+ 
into ([30]) we get 

7 / N + 7+ ^ + 7+ .X 

^ 53 + a52+5(/3(0)+/3++/3+_)+7++7+_ (8 + a;i)(5 + a;2)(5 + a;3) ^ ^ 

where cJi ,CJ2 and Us are roots of the following equation 

uj^ - auj^ + uj{P{0) + /3+ + /3+_) - (7+ + 7+_) = (35) 

Inverse Laplace transformation of equation (|34|) gives the exact expression of p++?/;++(t) 

e--^^(7+-/3+c^i) , e-^^\j^ - p^uj2) , e-"3t(^^ _ ^^^3) 
P++^++(^) = 7 ^7 V + 7 u V + 7 u V (^^) 
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Similarly, using the expressions of ao and a+ in eq. (jSTJ), we get 

- ^3 ^ ^ ^^^^^^ ^ ^ ^ ^ (d7j 

Inverse Laplace transformation gives the exact expression of 

I 
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FIG. 3: (Color online) Probability density of conditional dwell 
time, — (t) is plotted for a few different values of parameter 
(jjr. The values of other parameters are (all in s~^) C(;/=20.0, 
cjp/^SO.O, cjpr^lS.O, c<;^=40.0 and Ue =4.0. 



FIG. 4: (Color online) Probability density of conditional dwell 
time, '0++(t) is plotted for a few different values of parameter 
(jjf. The values of other parameters are (all in s~^) cJr^S.O, 
cjp/^SO.O, cjpr^lS.O, Uh^^O.O and Ue — 4.0. 



Note that by putting s = in equation (|34|) and ([37j) 
we get the "branching probabilities" 



7+ 



7+ + 7+- 



7+- 



7+ + 7+- 



(39) 



(40) 



which satisfy the normalization condition p++ +j9^ = 1. 

The cDTDs V^+_(t) and ^/^^^(t) are plotted in figs. [3] and 
m respectively, for a few different values of the parame- 
ters ujr and uof. In both the figures, the most probable 
dwell time increases with decreasing uOj. and decreasing 

UJf. 

In the same matrix-based formalism, the average ve- 
locity of a DNAP is given by the general expression [11 1 



.7(0) 



(41) 



where the dot indicates derivative with respect to q. The 
right hand side of eqn. (|4T]) can be evaluated for our model 
of DNAP using the explicit expressions ([19]) and (p!6|) for 
^{q) and /3(0), respectively. For the purpose of show- 
ing close relation of V with the Michaelis-Menten (MM) 



equation for the average rates of enzymatic reaction, we 
now assume that dNTP binding is rate limiting (the gen- 
eral framework of our theory does not need this assump- 
tion). Under this assumption, we can write 

ujf= u^fldNTPc] and 1]/ = n^f[dNTP^] (42) 

where [dNTPc] and [dNTP^j] are the concentrations of 
the correct and incorrect substrates, respectively, and 
that uj^ » Q^j. In this case, the average velocity of 
the DNAP, i.e., the average rate of polymerization, can 
be expressed in a MM-like form 



^ K^atidNTP,] 
^ Km + [dNTP, 

for the correct nucleotides, where 



Kcat — 



and the effective Michaelis constant is 



Km = 



(43) 



(44) 



(45) 
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Replacing uj by Q and [dNTPd by [dNTP^] we get the 
average rate of polymerization for the wrong nu- 
cleotides. In the limit of negligible exonuclease activity, 
the kinetic diagram shown in fig (T] reduces to the scheme 
shown on the left panel of figE] which is the standard 
MM-scheme with a single intermediate complex; in this 
limit the expressions for Kcat and Km are consistent with 
those for the standard MM scheme |l2| . 



B. Distribution of turnover time for exonuclease 

In this section we derive the DTT for exonuclease ac- 
tivity of the DNAP. We insert a hypothetical state 
such that 



where 



Ps ^ Pi* 4 Pi 



(46) 



where in the limit S ^ oo, Pi and Pf become identical 
and we recover our original model. 

For the simplicity of notation, in this subsection we drop 
the site index without loss of any information. The mas- 
ter equations for P/i(t) (/i = 1,2,3) and that for Pi{t) 
are 



dPijt) 
dt 



dP2{t) 
dt 



dPsjt) 
dt 



= -UJfPi{t)^{Ur^UJh)P2{t) (47) 

ujprPsit) (48) 



= UJpfP2{t) - {uJe + UJpr)P3{t) 



(49) 



N 



and 



—UJf UJh -\- (jOr 

UJpr 

UJpf -(We+Wpr) 



Q 



Piit) 

P2{t) 

Ps{t) 



(53) 



(54) 



Solution of the equation (|52]) in Laplace space 

Q{s) = S(5)-^Q(0) (55) 

where 

S{s) =sl-N. (56) 
Solution (|55|) for the assumed initial conditions provide, 

(-l)l+3Coi3 



Psis) = 



|S(s)| 



(57) 



which, explicitly in terms of the rate constants, takes the 
form 



p rsi = ^/^p/ 



(58) 



where 



a' = UJe -\- UJf -\- UJh -\- ^pf + ^pr + (59) 



dt 



.Psit) 



(50) 



For the calculation of DTT, we impose the initial con- 
dition Pi(0) = 1, and P2(0) = P3(0) = Pi*(0) = 0. 
Suppose, f{t) denotes the DTT. Then, f{t) At is the 
probability that one exonuxlease cycle is completed be- 
tween t and t + At, i.e., the DNAP was in state 3 at time 
t and made a transition to the state 1* between t and 
t + A. Obviously, fit) At = uJeP3{t) and, hence. 



f(t) = LO.Psit) 



(51) 



Using a compact matrix notation, the equations (|47|) , (|48 ]) 
and (|^9]) can be written the form 



-Q{t) = NQ(t) 



(52) 



/3' = UJf UJpf + UJeUJh + UJeUJpf + ^f^pf + UJf UJpr + UJhUJpr 
+ UJeUJpr + UJpruj^ (60) 

y=UJfUJpfUJe (61) 

Since in the Laplace space the equation (|5T]) becomes 

/(S) = LUePsis), (62) 

we get 

i-/ X ^e^/^p/ i^e^f^pf 

^ ^ s3 + a's'^ + 13' s + 7' (s + ^;l)(s + ^;2)(5 + v^) 

(63) 

as the DTT in the Laplace space. 
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Taking the inverse Laplace transform of equation we get the DTT 



/(*) 



UJfUJpfUJe 



_{vi - V2){vi - Vs) 



UfUpfUJe 



{V2 - Vi){v2 - Vs) 



{vs - vi){vs - V2) 



-vst 



(64) 



where vi^V2,vs are solution of the fohowing equation 



{uJe ^ UJf ^ UJh ^ UJpf + UJpr + UJr)v'^ + {^f^pf + ^e^h + ^e^pf + ^f^pf + CJfUJpr + UJhUJpr + CJg^pr + ^pru;,)^ 



(uJfUJpfUJe) = 0. 



(65) 
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FIG. 5: (Color online) Probability density of turnover time 
for unproductive exonuclease mode of enzyme, f(t) is plotted 
for a few different values of parameter ujpf . The values of the 



other parameters are (all in s ^) C(;/=20.0, Ur 
cj^=40.0 and uje — 4.0. 



=5.0, C(;pr=15.0, 



FIG. 6: (Color online) $p and $e plotted against ojpr while 
the ratio Qpr/oopr — 0.1 is kept fixed. The three curves cor- 
respond to cje = 1.0s~^, 2.0s~^, 3.0s~^. The values of the 
other pararemeters are (all in s~"^): ujf — 1.0, = 10~^, 

UJpf 0.1, VLpf 10.0, UJr .1, VLr 1.0, UJh 10.0, 

r^/, = 0.1. 



This DTT is plotted in fig(j5j). Plots are consistent with 
the intuitive expectation that increasing ujpf leads to a 
decrease the turnover time. 

Suppose {t) denotes the mean time gap between the 
completition of the successive exonuclease reactions cat- 
alyzed by the DNAP. Then the average rate = 1/ (^) of 
the exonuclease reaction can be expressed in a MM-like 
form 



K,at[dNTP,] 



with Knat 



UJpfUJe 



Km + [dNTPc\ " UJe^UJpf ^ UJpr 

(66) 

for the correct nucleotides where Km = Km- Replacing 
uj by n and [dNTPd by [dNTP^] in (j66j) we get Ve^"^^ 
for the wrong nucleotides. In the limit uJh ^ 0^ ujpr 0, 
the kinetic diagram shown in fig (1] reduces to the simpler 
scheme shown on the right panel of figl2] which is essen- 
tially a generalized MM-like scheme with two interme- 
diate states. Not surprisingly, in this limit, the average 
rate of the exonuclease reaction is consistent with that of 
the MM-like scheme with two intermediate states [12i] . 



C. Quantitative measures of error 



Note that ^p = Vp^'^^/iVp^'^^ 



Vp^"^) is the fraction of 



nucleotides misincorporated in the final product of repli- 
cation. Similarly, the fraction = / i^e + K^) ^ 
measure of the errorneous severings^ i.e., fraction of the 
cleaved nucleotides that were incorporated correctly into 
the growing DNA. Since ujpr is the strength of the "cou- 
pling" between the two different enzymatic activities, we 
plot ^p and against ujpr in figl6]for a few typical sets 
of values of the model parameters. 

Decreasing $e with increasing ujpr is a consequence 
of the escape route via uOpr for the correctly incorpo- 
rated nucleotides that get transferred unnecessarily to 
the exonuclease site. It is the increasing number of such 
correctly incorporated nucleotides recused from the ex- 
onuclease site that leads to the lowering of ^p with in- 
creasing UJpr. The limiting values of ^p and in the 
limit of large UJpf , are determined by the corresponding 

limiting expressions Vp^-' {ujfUJh)/{uJh +^/ +^r), and 

Ve^-' 0:^ {UJ f UJpfUJe) /[uJpr{0Jh +^r)]- (ExprCSSioUS for 
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vr^ and vr^ are similar in the limit of large Vtpr-) 

IV. SUMMARY AND CONCLUSION 

Here we have theoretically investigated the effects of 
the coupling of two different modes of enzymatic activ- 
ities of a DNAP, in one of these it elongates a DNA 
whereas in the other it shortens the same DNA. The fun- 
damental questions we have addressed here in the context 
of DNA replication have not been addressed by earlier 
theoretical works [14]. The effects of tension on the poly- 
merase and exo nuclease activities, which have been the 
main focus of the earlier works [14], will be reported else- 
where [l^ . The mechanism of error correction by DNAP 
is somewhat different from the mechanism of transcrip- 
tional proofreading which is intimately coupled to "back 
tracking" of the RNA polymerase [16, Tz]. 

We have derived exact analytical formulae for the 
cDTD and DTT which will be very useful in analyzing 
experimental data in single DNAP biophysics, particu- 
larly its stepping patterns and enzymatic turnover. In 



spite of their coupling, the average rates of both the 
enzymatic activities are MM-like; the analytical expres- 
sions for the effective MM parameters explicitly display 
the nature of the coupling of the two kinetic processes. 
We have also reported exact analytical expressions for 
the fractions and which measure replication error 
and erroneous cleavage] these expressions can be used 
for analyzing data from both single molecule [5] and 
bulk [18| experiments on wild type and mutant DNAPs. 
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